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Abstract 

We derive the high energy asymptotic of one- and two-loop corrections in the next- 
to-next-to-leading logarithmic approximation to the differential cross section of H^-pair 
production at the LHC. For large invariant mass of the W-pair the (negative) one-loop 
terms can reach more than 40%, which are partially compensated by the (positive) 
two- loop terms of up to 10%. 
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1 Introduction 



With the LHC now starting its operation, the experimental investigation of scattering pro- 
cesses at the TeV scale is within reach. Starting from these energies electroweak corrections are 
strongly enhanced by Sudakov logarithms of the form (a"^ In^" s/M^)" [Ij. The full evaluation 
of electroweak one-loop corrections to fermion- or VF-pair production is by now a straightforward 
task. Two-loop corrections, however, can be obtained only in the high energy limit. By employing 
the evolution equation approach the analysis of the dominant logarithmically enhanced two-loop 
corrections for four-fermion processes has been pushed successfully from next-to-leading logarith- 
mic (NLL) approximation [2|3I4] to NNLL [5] and even N^LL approximation [6|7j . which accounts 
for all the two-loop logarithmic terms (for additional work on this topic see e.g. }8|9|10] ). Subse- 
quent analysis performed in the effective theory framework [llj employing the two-loop anomalous 
dimensions calculated in Refs. |6|7] have confirmed the formentioned result. 

In this paper we consider specifically pair production of VT-bosons. Previously, the electroweak 
corrections were studied mainly in the context of the electron-positron annihilation. The one-loop 
corrections have been evaluated for the W^-pair production [12|13|14|15j and the VF-boson mediated 
e^e~ — )• 4/ processes jl6ll7ll8ll9j . For high energies the two- loop logarithmically enhanced terms 
have been obtained up to the NNLL approximation |3I4|20|21] . The one-loop contribution amounts 
to typically -20% for 1 TeV and -50% for 3 TeV while two-loop terms vary between 2 and 5% 
for 1 TeV, for 3 TeV they may even rise to 20%. For the I^-pair production at the LHC the 
analysis of the one-loop electroweak logarithms to the NLL approximation is given in |22|23j with 
the realistic cuts and the effect of gauge boson decay included. Beyond one loop the logarithmic 
corrections to the partonic cross sections were considered in [24]. In view of the extremely large 
partonic energies and with the LHC eventually operating at full luminosity (not to speak of the 
SLHC) invariant masses of the Il^-pair exceeding 1 TeV and approaching 3 TeV seem within reach. 
Therefore the evaluation of the enhanced electroweak corrections is of particular interest. Here we 
present the explicit result for the one- and two- loop corrections to the partonic qq — )• W~^W~ and 
hadronic pp — )• W~^W~ cross section in high energy limit in the NNLL approximation. Note that 
the cross section of W^-pair hadronic production is a subject of large corrections due to the strong 
interaction of the initial states. Currently the analysis of QCD corrections is completed to the NLO 
and NLL approximation (see [25|26|27] and references therein). The size of the corrections depends 
strongly on a particular observable and in many cases the available approximation provides a few 
percent accuracy. As we will see the two-loop electroweak logarithms become essential at this level 
of precision and have to be included in the theoretical predictions. 

Our paper is organized as follows: the partonic processes in Born approximation are introduced 
in Section [21 In Section [3] the evolution equation approach is outlined for the simplified case of a pure 
SU{2) spontaneously broken gauge theory. The discussion closely follows Ref. [21]. However in the 
present paper we derive the explicit result for the one-loop corrections to scattering amplitudes given 
in Appendix A. The generalization to the SU{2) (8) U{1) Standard Model is presented in Sectional 
which contains a more detailed analysis of the separation of infrared singularities connected with 
virtual photon emission. The results for the one- and the two-loop corrections to the partonic 
cross section in NNLL approximation are listed in the Appendix B. In Section 14.21 we present a 
numerical study of these corrections for ^/§ = 1 TeV and 3 TeV respectively. Based on these results, 
the corrections to the transverse and longitudinal Vl^-pair production in proton-proton collisions at 
14 TeV are presented in Section [431 together with the discussion of the anticipated statistical errors. 
Section [5] contains a brief summary and conclusions. In Appendix C we present the correction to 
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the two-loop NNLL result for the transverse T^-pair production in electron-positron annihilation 

2 The partonic process 

The partonic processes relevant for the VF-pair production at hadron colliders are gluon fusion 
and quark-antiquark annihilation. The gluon contribution to the total cross section is about 5% 
[28] and we focus on the process qq ^ W~ . In the leading order it is described by the diagrams 
in Fig. [TJ 




Figure 1: Tree level diagrams contributing to the partonic process 
The kinematics at partonic level is defined by: 

q{pi, A+) + q{p2, A_) ^ W+{k+, K+) + W~{k^,K-), (1) 

where X± and k± are the helicities of the incoming and outgoing particles respectively. For on-shell 
VF-bosons, the matrix element can be then expressed as function of the Mandelstam variables: 

s = {pi+P2f, i={pi-k^f, u={pi-k+f, (2) 

They are related to the scattering angle 9 through the relations: 

i = M^--^(l-P cose), n = M2,-|(l + /3cos0), = i _ 4 ^ . (3) 

In the high energy limit only final states where the VF-bosons have the same polarization are not 
suppressed by a factor M^/s or higher. In addition, the case where both VK's are longitudinal can 
be reduced by means of the Goldstone equivalence theorem to the production of a pair of charged 
Goldstone bosons as shown in Fig. [2j 




Figure 2: Goldstone equivalence theorem at Born level 

3 Massive gauge boson production in SU{2) model 

Let us, in a first step, neglect the hypercharge and consider a simplified model with sponta- 
neously broken gauge group SU{2). The model retains the main features of the massive gauge 
boson sector of the Standard Model. In this case the result can be presented in a simple analytical 
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Figure 3: Fermion/scalar scattering in an external singlet vector field and scattering of a gauge boson in an 



external scalar field. The momentum of external field satisfies -p 



form and constitutes the basis for the further extension to the full electroweak theory. We study 
the process of gauge boson pair production in fermion-antifermion annihilation at high energy and 
fixed angle with all kinematical invariants of the same order and far larger than the gauge boson 
mass M, |s| ~ \t\ ~ \u\ ^ M^. In this limit the asymptotic energy dependence of the amplitudes 
is dominated by Sudakov logarithms |29|30] and governed by the evolution equations [3 1|32|33] . 
The method of the evolution equations in the context of the electroweak corrections is described in 
detail for fermion pair production in Ref. [7] and for PV^-pair production in Ref. [21] . 

Following Ref. [21j we introduce the functions Z^,(j,,A which describe the asymptotic dependence 
on the large momentum transfer Q of the scattering amplitude of the spinor (ip) or scalar {(p) field 
in an external singlet vector field and of the vector boson (^4) in an external singlet scalar field, 
i.e. of the respective form factors in the Euclidean region (see Fig. [3]). In leading order in M^/Q^ 
these functions are known to satisfy the following linear evolution equation |31|32|33] 



d 



51nQ2 



' —j,{a{x)) + aa{Q^)) + U»{M^)) 



with the solution 




exp< / — 



— 7.(a(x')) + C^(«(^))+6(«(M2)) 

M2 X' 



(4) 



(5) 



which satisfies the initial condition -Zi|Q2_jvf2 = 1- Here the perturbative functions 7i(a) etc. are 

given by the series in the coupling constant a(/x^), e.g. 7i(a) = Y^^=i{^/^'^)'^li^^ ■ Then the 
amplitude of the transverse (longitudinal) gauge boson production At {-Al) can be decomposed 
as follows 

At,l = ol{^t,l)2,^Za,cI,At,l ■ (6) 

where A.t,l is the reduced amplitude and we factor out the Born coupling constant a{fj,T,L)- The 
scale dependence of this factor is cancelled by the higher order renormalization group logarithms 
replacing fiT,L by a physical scale of the process. In the case of the longitudinal VF-pair the proper 
scale is fiL = ^/s because it describes the interaction of far off-shell intermediate gauge boson with 
virtuality of order ^/s. For the transverse PV^-pair production it is /u^ = Mw corresponding to the 
coupling of the on-shell VF-bosons. Note that in an alternative approach based on the soft-collinear 
effective theory [24] the normalization scale of the Born coupling constant for the transverse gauge 
bosons is set to ^/s. This is compensated by an additional /3o contribution to the anomalous 
dimension (^^'^ which effectively shift the normalization of the Born coupling constant to M^y, in 
agreement with our result. 

Due to the factorization property of the Sudakov logarithms associated with the coUinear diver- 
gences of the massless theory [35j the reduced amplitude satisfies the simple renormalization group 
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like equation [36|37|38] 



jA.T,L = XT l("(<3^))^T,L 



91nQ2'-'- -^J,^^-^^ (7) 

where = —s and Xt L is the soft anomalous dimension matrix acting in the space of the isospin 
amplitudes. The solution of the above equation is given by the path-ordered exponent 



A 



■T,L 



Pexp 



Q da; 

—Xt,l{ol{x)) 

M2 X 



A)T,l(«(M2)), 



(8) 



where Aqtl determines the initial conditions for the evolution equation at Q = M. By calcu- 
lating the functions entering the evolution equations order by order in a one gets the logarithmic 
approximations for the amplitude. By expanding the exponents one gets the one- and two-loop 
corrections in the following form 



(1) 



A^-J, + a!^tt, 



+ - 
2 



-/3o(C«+xgl 



7(^) = 7?^+7il, (9) 



where L = \t\.{Q'^ / M"^) and /3o is the one-loop beta function. The anomalous dimensions 7(a), 
C{oi) and x(a) are mass- independent and can be associated with the infrared divergences of the 
massless (unbroken) theory. At the same time the functions £,i{oi) and Aq^t,l{o) do depend on the 
infrared structure of the model and require the calculation in the spontaneously broken phase. All 
the perturbative coefficients in Eqs. ([9]) except .a!^\, ^ are known [21]. In Ref. [21] the result for the 
one-loop correction to the cross section [19] has been used to obtain the two-loop NNLL terms. We 
complete this part of the calculation and present the explicit result for the one-loop corrections to 
the amplitude in Appendix A. Our result for the cross section agrees with Ref. [19j . 

The large Yukawa coupling of the third generation quarks to the scalar (Higgs and Goldstone) 
bosons results in specific logarithmic corrections proportional to m^/M^. This kind of Sudakov 
logarithms were studied in Ref. [32j and have universal structure for any renormalizable non-gauge 
theory. The factorization in this case is much simpler than in gauge theories and the logarithmic 
corrections are completely determined by the ultraviolet field renormalization of the external on- 
shell lines. Since the Yukawa coupling of the initial light quark states is suppressed the Yukawa 
enhanced Sudakov logarithms for hadronic production of VF-pair are similar to those for VF-pair 
production in electron-positron annihilation [21]. Thus the Yukawa enhanced corrections can be 
taken into account through the modification of the evolution equations for the corresponding 2^^- 
function. The main complication is that the Yukawa interaction mixes the evolution of the quark 
and scalar boson form factors and in general does not commute with the SU{2) coupling. Thus 
the evolution equation has a complicated matrix form: 



91nQ2 



dx 



7(a(x)) + C(a(Q'), ay«fe(Q')) + i{(^{M^)) 



(10) 
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with the solution 




^ dx' 

— 7(a(x')) + C(a(x), Qy„fe(x)) + $(,a{M^)) 
A/2 x' 



(11) 



where 7^^) = (-3/2)1, ^ = 0, ayuk = / (2M^) a, and we introduce the five-component vector 

Z = {Z^, Z^, Zb^, Zt^, Zt+) . (12) 

The subscript + (— ) stand for the right (left) quark fields and Z^ corresponds to the transition of 
the Higgs boson into the neutral Goldstone boson in the external singlet vector field. The one-loop 
anomalous dimension matrix reads [2T| 



/12 0\ 
12 
9 
9 

VO 0/ 



/ 6 -6\ 



+ 





1 
1 

V-1-1 



6-6 
0-1 
0-1 
-1-1 / 



(13) 



where the first term represents the pure 5'C/l(2) contribution, the second term represents the 
Yukawa contribution and we introduce the ratio p = ayuk/oi = .A ~ 1- The proper initial 



condition for the evolution equation which corresponds to the SUl{2) Born amplitudes of the third 
generation quark and scalar boson production in light quark-antiquark annihilation is given by the 



vector Z 







2(r|, 



T^3 



1^1,— 1,-1,1,0) where T? stands for the particle isospin 



and the overall factor of 2 is introduced for convenience. Since the Yukawa enhanced logarithmic 
corrections can be attributed to the external on-shell field renormalization we expect a diagonal 
form of the corrections. This is indeed the case due to a nontrivial matrix relation 



(1) f " 

Yuk I 



■Zn 



2 



^Yuk J 



2n+l 



■Zn 



3p2 
2 



^ {-3p,3p,p/2,-p/2,0) , (14) 



where the vector on the right hand side of the second equation represents the one loop correction 
Cylk ' ^0- By factorizing the components of Zq we can rewrite it as follows 



CiL-^o = 2pi-3Tl-3Tl 



1. 



1, 



-T; 



(15) 



The coefficients of in the above expression depend only on the field renormalization of the particle 
i as it has been explicitly shown in Ref. [39j. For example, for the hypercharge mediated Born 



amplitudes of the same process we have different initial conditions Zq = {Y^, Y^, yj,„, Yi_, Yt 



t+) 



(1, 1, 1/3, 1/3,4/3) where Yi stands for the particle hypercharge. However the one loop coefficients 
are the same as in Eq. (jl5p 



^Yuk ' ^0 



p(-3y^,-3K„-in_ 



1 



Yt-,-Yt 



t+) 



(16) 



The different form of the odd and even order corrections is dictated by the off-diagonal character 
of the matrix of field renormalization by Yukawa interaction. 
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By expanding the solution for the component we obtain the Yukawa enhanced logarithmic 
corrections to the amplitude of the longitudinal 1/F-pair production. Let us introduce the following 
notation 

{Ovuk = [C • Zo]^ , (17) 

where only the terms proportional to the second or fourth power of the top quark mass are kept on 
the right hand side. Then the Yukawa contribution to the amplitude ([9]) takes the following form 



Yuk 



^ I Yuk 



ly'HC^'^)YukL' + 



(1)1 



Yuk 



Yuk, 



Yuk 



(18) 



where /3|^"'^ = 9/4—3/9/2 is the one-loop beta-function of the Yukawa coupling constant and a!^^ {vuk 
is the one-loop nonlogarithmic Yukawa contribution given in the Appendix A. 



4 W-pair production in the electroweak model 

4.1 Analytic results 

The electroweak Standard Model with the spontaneously broken SUl{2) x U{1) gauge group 
involves both the massive W and Z-bosons and the massless photon. The corrections to the fully 
exclusive cross sections due to the virtual photon exchange are infrared divergent and should be 
combined with real photon emission to obtain infrared finite physical observables. The infrared 
divergences of the virtual corrections are regulated by giving the photon a small mass A. In the 
limit <C My^r <C the dependence of the amplitudes on A in the full theory is the same 
as in QED. Thus the logarithmic corrections can be separated into "pure electroweak" Sudakov 
logarithms and QED Sudakov logarithms of the form ln(Q'^/A'^) or ln(M^/A^). 

To disentangle the electroweak and QED logarithms we use the approach of Ref. |1|5|7] . While 
the dependence of the amplitudes on the large momentum transfer is governed by the hard evolution 
equations (c./. Eqs. (jHE])), their dependence on the photon mass is governed by the infrared evo- 
lution equations flj . Two sets of equations completely determine the dependence of the amplitudes 
on two dimensionless variables Q/Mw and Q/X up to the initial conditions which are fixed through 
the matching to the fixed-order result. For A^ <^ Myy the singular dependence of the amplitudes on 
the infrared regulator is governed by the QED evolution equation. Its solution to NNLL accuracy 
in the massless fermion approximation = {f ^ t) is given by the factor 



U= Uo{ae)exp 



(47r)2 9 



+ 1) ^ + f 3Ql - 4Q, InJ ^ + In' ^ + 2 In 



ae(A" 
47r 



" ' ; A2 ' V -^'?-£y-A2 ' A2 ' A2 

+ 0{4)\, (19) 



X2 ' 3 i A2 ' V A2 



where is the MS QED coupling constant, and Qq is the quark electric charge. The NNLL 
approximation for U can be obtained from the result for the fermion- antifermion production [5J by 
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proper modification of the QED anomalous dimensions. Note that we take into account the top 
quark decouphng and Eq. (jl9p corresponds to five Ught flavors in contrast to Ref. [5] where all the 
quarks were assumed to be massless. To exclude the top quark contribution in the expressions for 
the QED anomalous dimensions in [5] Ng should be replaced by Ng — 1/2. The preexponential 
factor Uq in Eq. (I19p is factorization scheme dependent. It is convenient to fix it by normalizing 
^(ae)| ,^x2=M2 ~ ^' factorize the QED factor and write the full theory amplitude as a product 



A= UAe 



(20) 



where includes only electroweak Sudakov logarithms. The logarithms of the photon mass in 
U are generated by loops with soft photons, photons collinear to the initial state fermions, and 
soft photons collinear to the final state gauge bosons, which result in the logarithmic dependence 
of the coefficients on My/. In the physically motivated cross section which is inclusive in respect 
to the photons with the energy much less than electroweak scale the singular dependence of lA 
on the photon mass is replaced by the experimental cuts on the soft photon energy or absorbed 
into the parton distribution functions. One may easily change the regularization scheme and use 
e.g. dimensional regularization which is more convenient for the analysis of the parton distribution 
functions. In the present paper we focus on the pure electroweak part of the amplitude Aew The 
factorization formula (I20p implies that the anomalous dimensions corresponding to the electroweak 
Sudakov logarithms are obtained by subtracting the QED contribution from anomalous dimensions 
of the full theory. The functions 7, and % are mass- independent. Therefore the anomalous 
dimensions parametrizing the electroweak logarithms can be obtained by subtracting the QED 
contribution from the result of the unbroken symmetry phase calculation to all orders in the coupling 
constants. In particular in one loop we get 



'A,(t> 
^(1) 



X? 



SU(2) 



SU(2) 



+ 4g,4ln^l. 



^(1) 



- ^A,(l} 



SU{2) 



SU(2) 



V2 



Y,Y.tl + AQ,s^ 



1 '^w 



u 

ln-1, 21 

t ^ ' 



where K, Ya = 0, 1^ 



-1 are the hypercharges of quarks, gauge and Goldstone bosons, 



sin 9y/, = tan 9w, and 9w is the electroweak mixing angle. The SU{2) part of the anomalous 
dimensions can be found in [21] while the hypercharge contribution and QED subtraction term are 
given explicitly. The anomalous dimensions for the quark 2^- functions can be found in Ref. [5] . The 
only two-loop coefficients we need are 



(2) _ (2) 
'^A,</. ~ ^A4 



SU(2) 



800 



— QA,<j,S^ : 



(22) 



in the MS scheme. On the other hand the functions and ^0 are infrared sensitive and require 
the use of the true mass eigenstates of the Standard Model in the perturbative calculation. In 
NNLL approximation one needs the one-loop contribution to these quantities which can be found by 
comparing the solution of the evolution equation with the explicit one-loop result for the amplitudes. 



In this way we find that the anomalous dimensions ^, 
difference between M]y and Mz and obtain: 



(1) 



get contributions just from the mass 



(1) 



■3\2 



2) ^■ 



In 



(23) 
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where Tf = Qj — 1^/2 is the third component of the isospin. The expressions for the nonlogarithmic 
one-loop corrections to the amphtude ^q^^ are rather cumbersome and we collect them in Appendix 
A. Note that ^q^"* depends on the normalization of the QED factor. We use the normalization where 
all the nonlogarithmic one-loop corrections are contained in .4q^''. With the above parameters of 
the evolution at hand we can write down the two-loop NNLL corrections to the amplitudes as in 
Eq. Q. The two-loop Yukawa contribution in the NLL approximation is given by the interference 
of the one- loop double logarithms and the one-loop Yukawa enhanced single logarithms. Thus it is 
straightforward to obtain this contribution exactly. For the NNLL two-loop Yukawa contribution 
we use the SU{2) model of the previous section with p = m^/(2M^), which approximates the 
exact result with the accuracy of order sin^ Oi^ ~ 0.2. 

Now we are in the position to present the final result for the cross sections. We define the 
perturbative series as follows 



The coefficients for the one and two-loop NNLL terms are listed in the Appendix B. Below we 
present the numerical analysis of the corrections to the partonic and hadronic cross sections. 

4.2 The partonic cross section 

For the numerical estimates we adopt the following input values 



Mw = 80.41 GeV, Mz = 9L19 GeV, Mh = 117 GeV, mt = 172.7 GeV, (25) 



and take = 1 TeV as characteristic example. The one and two-loop corrections for left-handed 
li-quarks in the initial state are plotted in Fig. U] showing a sizable NNLL contributioE0. The 
structure of the corrections for the left-handed d-quarks is similar, see Fig. [3 To facilitate the 
comparison of the u- and d-quarks cases related by crossing symmetry in the Born approximation, 
we plot the cross section for u-quarks as a function of — cos 9. In the Born cross section we always 
use the physically motivated normalizarion scale of the coupling constants, which is ^ = M]^ for 
the transverse and // = ^/s for the longitudinal boson production. 

The contribution of the right-handed quarks vanishes for transversally polarized T^-bosons, and 
for longitudinally polarized bosons it is significantly smaller than the one of left-handed quarks, see 
Fig. [6l In one as well as in two-loop approximation one observes large compensations between LL, 
NLL and NNLL terms. Evidently the LL approximation, even when combined with NLL terms 
only, does not lead to an adequate description of the full result. In Ref. [T5] the quality of the high 
energy approximation has been studied at one-loop level. The error turns out to be less than a few 
percents for a partonic center of mass energy above 500 GeV and a scattering angle in the range 
30" < 6* < 150°. 

^ Numerical results for the partonic cross section have been presented in Ref. [24) and qualitatively agree with 
our analysis. However a direct comparison of the results is not possible since the authors of [24] use different power 
counting and QED subtraction prescription. 




(24) 



a(M|) 



1 



si, = 0.231, 



128.1 
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Figure 4: One and two-loop corrections to the partonic cross section for left-handed w-quarks in the initial 
state, transverse (left panel) and longitudinal (right panel) M^-bosons at = 1 TeV. 



4.3 Hadronic cross section 

To obtain transverse momentum and invariant mass distributions for the process pp — t- W^W~ + 
X the partonic cross section must be convoluted with the parton distribution functions fhi,i{xi,lJ'p) 
and fh2,jix2, fi"^), where fj,F is the factorization scale, xi and X2 are the momentum fractions 
carried by the parton i in the hadron hi and by the parton j in the hadron h2 respectively. The 
pr-distribution is given by 



da 



1 d& 

jV2 Jf^^^^ J^dX2fhi,iixi,l4^) fh2j{x2,IJ,F) 9{xiX2- Train) (26) 



dp 

where Nc is the number of colors, the sum is over all possible qq partonic initial state, pr = 
sm9^s — 4:Myy 1 2 is the transverse momentum of the T^-bosons and we adopt [ip = Pt- The 
quantity 

s 

is related to the minimal partonic energy that is needed to produce two VF-bosons with a given 
transverse momentum px- The partonic differential cross section d&ij/dpT are given in terms of 
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Figure 5: One and two- loop corrections to the partonic cross section for left-handed d-quarks in the initial 
state, transverse (left panel) and longitudinal (right panel) Vt^-bosons at Vl — 1 TeV. 



the angular differential cross section as follows 
d&ij 4pT 



S Try 



-f^^ + {t^u) , s = X1X2S. (28) 
a cos a 



The numerical results are obtained by using the MRST parton distributions [30] and the integra- 
tion routine CUHRE from the CUBA library [51]. The upper panel of Fig. [7] shows the NNLO 
PT-distributions for the production of transverse and longitudinal VF-bosons in the NNLL approx- 
imation. Transverse bosons production is evidently dominant, with the cross section being about 
twenty times larger than the one of the longitudinal bosons. The lower panel of Fig. [7] shows the 
NLO and NNLO corrections separately. For the production of transversely polarized VF-pairs the 
one-loop correction reaches 40% at px = 1 TeV and 60% at = 2 TeV. The two-loop contribution 
amounts up to 10% at px = 1 TeV and 20% at px = 2 TeV and partially compensate the one-loop 
corrections. For the longitudinal boson production the one-loop correction is about 15% (30%) at 
px =1 TeV (px =2 TeV), while the two-loop contribution does not exceed a few percent up to 
px = 2 TeV. As anticipated above the radiative corrections for the longitudinal case are smaller 
than those for transverse W bosons. This is because the value of the quadratic Casimir operator 
of the SUl{2) electroweak group, which govern the leading logarithmic contribution, is smaller 
for the fundamental representation of the longitudinal degrees of freedom than for the adjoint 
representation of the transversely polarized M^-bosons. 
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Figure 6: One and two-loop corrections to the partonic cross section for right-handed u-quarks (left panel) 
and d-quarks (right panel) in the initial state, and longitudinal M^-bosons at \fh = 1 TeV. 



The invariant mass distribution for the PF-pair production is defined as follows 



da 



-y 

I] 



^ f ( 2^f f 2Nd'^ii(M^ a) 

dX2jhi,i{^i^f^F)Jh2,j{X2,f^F)- 



dM„ 



(29) 



where Mj^^ = (A;+ -|- A;_)^ = \/I is the invariant mass of the VF-pair system and we adopt 
= M^^. Here the partonic differential cross section daij/dMy^^^^r is obtained by integrating the 
angular differential cross section in the region — cos 6'min < cos 6 < cos 6'min 



daij{M' 



2 •■ 
WW ' 



dM„ 



d cos Q- 



d<Ti,-(M 



2 ^ 

WW > 



dcos0 



(30) 



which excludes the range of small angles where the high energy and the Sudakov approximations 
are not valid. The results for the invariant mass distribution are plotted in Fig. [8] with an angular 
cutoff ^min = 30°. To estimate the potential statistical sensitivity, the corresponding plots are 
shown for the production cross section of W pairs with px > p^*. Taking, as crude estimate, an 
integrated luminosity of 200 fb~^, about 1200 W pairs with pT > 600 GeV would be produced. 
Assuming that the experimental analysis would be based on the final state with one VF-boson 
decaying leptonically and the other hadronically a fraction of about 4/9 of the pairs could be 
observed, leading to a nominal statistical error of about 4%. Under this (optimistic) assumption 
the one-loop terms would be clearly relevant and the two-loop terms start to contribute. 
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Figure 7: Transverse momentum distribution (including corrections) of transverse and longitudinal W pairs 
and relative corrections for proton-proton collisions at ^/s — 14 TeV. 

5 Summary 

In the present paper we derived the one and two-loop electroweak corrections to T^-pair produc- 
tion at the LHC in NNLL approximation in high energy limit. We present the analytical result for 
the amplitudes, differential partonic cross sections, hadronic pT- and invariant mass distributions. 
The structure of the corrections is similar to the VF-pair production in e^e~ annihilation [21]. In 
the case of the transverse boson production we observe the cancellation between the huge NLL and 
NNLL contributions so that the sum is dominated by the LL term. For the longitudinal bosons the 
corrections exhibit significant cancellation between the LL, NLL and NNLL terms. The maximal 
effect of the corrections is on the pr-distribution of the transverse VF-pair production and reaches 
60% and 20% at px = 2 TeV in one and two loops, respectively. To push the theoretical error below 
1% the evaluation of the two-loop linear logarithmic terms should be completed, which requires the 
calculation of the two-loop mass-dependent anomalous dimensions 
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Figure 8: Total invariant mass distribution (including corrections) of transverse and longitudinal W^-pairs 
(upper panel) and the corrections to the invariant mass distribution (lower panel) for proton-proton collisions 
at ^/i = 14 TeV. 

Appendix A 

In this appendix we give the explicit analytical result, valid in the high energy limit, for the 
one and two-loop NNLL amplitudes of the processes: 

ATq_ : q^{pi) + q^{p2) ^ W}{k+) + VF^ (fc„) , 

: qTiPi) + <l±iP2) ^ + <P~ik~) , (31) 

where are the left/right handed fermions in the initial state, which can be either u or d quarks. 
The amplitude At describes the production of transversely polarized VF-bosons and vanishes for 
the right-handed initial quarks. The amplitude ^^i^ describes through the one-loop Goldstone 
equivalence theorem the production of the longitudinally polarized VF-bosons. The one-loop cor- 
rections to the Goldstone equivalence theorem of Fig. [2] can be properly described by introducing 
the following effective wave function counterterm for the 4>^ field: 

where is the longitudinal part of the W self-energy, S^*^ is the W—cj) self-energy, dM^r and 
6Z\Y are the mass and wave function counterterms of the W boson (see also 
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Figure 9: Cross section (including corrections) for transverse and longitudinal W pair production and 
relative corrections for proton-proton collisions at ^/s = 14 TeV (see text). 



The results given in this section are obtained by adopting the MS renormalization for the 
couphngs and the weak mixing angle and on-shell renormalization for the masses. As before the 
renormalization scale in the Born amplitudes is fixed to fi^ = for the transverse and ^\ = s 
for the longitudinal case. The Lorentz-Dirac structure of the amplitudes in the high energy limit 
takes a simple form: 



where uj± 



^Lq^ = q{p2) ^+ UJt liPl) -^Lq^ ' 



(33) 



-J— and we use the relation between the polarization vectors eK{k±) of the transversely 
polarized in the center of mass frame 



where k = it 1 stands for the polarization. The perturbative series for the amplitudes read 



n=0 ^ ^ n=0 



(34) 



(35) 
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where a{n^) = aeijji^) / s^{ii^) and the coupHng constants are supposed to be normahzed at = Mw 
unless the normahzation point is indicated exphcitly. The Born amphtudes read: 



(0) 



1 

U 



(0) 
Td- 



B. 



(0) 



0, 



(0) 



1 

2i 



1 + 



A 



(0) 



2s 



1 + 



^(0) ^ ,2 / 2 N 



(0) 



(36) 

The one-loop contribution to the second transverse Lorentz-Dirac structure of Eq. p3p is particularly 
simple and does not contain Sudakov logarithms 



B. 



(1) 



/-^T,S17{2) 



1 



(1) 



(37) 



i"t,SI/{2) 



9 n\ 5 
--5^ IL„, + -, 



1 / 1 n\ U r 9 nn 



1/3 u\ _^ 
9 V 2 t i 18 ■ 



(38) 



All the notations are explained at the end of the section. The one-loop corrections to the remaining 
Lorentz-Dirac structure can be formally decomposed according to the gauge coupling constant 
factor 



i(l) +/2 .(1) ,,2(^V 

•^Pq^, SU{2) ^ V •^Pq^^,Y ^ \-^Pq^, QED 



A 



(39) 



where the last term -^pq^ g^b corresponds to the first order term of the expansion of the singular QED 
factor (jl9p and cancels the infrared logarithms coming from soft photons and from photons collinear 
to the incoming quark^. The QED correction factorizes with respect to the Born amplitude and 
reads 



A 



(1) 

Pq^,QED 



2 ( Oj, + 1 ) L^fL + AQqL^Lut — QqL^ — 2 Q^L^L, 



'^q^l-^Z 



A 



(1) 



9 7 2 2 , r 

2~2^^^3^^ " 
9 



4 

A 



3 13 2 2,^ 



3q2+2)l^ + a(;]; 

L 



A 



(0) 



A 



(1) 

pgip ,sub 
2 . 



g^ + l)ln^^ + (3Q^-4Q,L„t 



A2 



In^ + ln 



2 ^ 2 In 



(40) 

where = —s — iO^ and Qq is the electric charge of the quark. Note that Jipq^^Q^r, vanishes for 
Mz ~ A — )• 0. After the subtraction the A-dependence disappears and we get the QED contribution 
in terms of the parameters of the evolution equation 



■^Pq^,QED -^pq^^sub 



(1) 

n ^ ^ \ Sq,QED 



,^(1) , ^(1) 
Zq,QED A 



.(1) 



q,QED I ^ + ^Pq,QED 



^In contrast to Section [4.11 we normalize here W(ae)j^2_^2_jv#2 = 1- In this case the amphtude is manifestly A 
independent. The numerical estimates are obtained with the normalization of Section [4. II where the A dependence 
survives in the imaginary part of the one-loop amplitude, but does not contribute to the cross section up to N'^LL 
approximation. 
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'lq,QED — + Cq.QBD — 3(5g, ^i,QED ~ ^Qq^ut^ Cg,QBD — 2 + 1 j -C/^ , 

^%QED = -[^QI+2)l, + AQ^L.Lut - Q\Ll - (41) 

In order to present the result of the remaining SU{2) and Y components, in terms of the coefficients 
of the evolution equations, it is necessary to analyze the isospin structure of the amplitude for left 
handed quarks. The general SU{2) basis for the amplitude Aj- of the left-handed quark-antiquark 
pair transition into two transverse gauge bosons reads 

q.q.^B^B,: (n+ d+) (^^i + ^2 + ^3 ^a^) (^J^ ) , a,6=l,2,3, (42) 

where Ba are SU{2) gauge fields and aa are the Pauli matrices. From the definition = 
{Bi =F iB2)/V2 we get the following structure for the production of W^W^ 



Ql Ql ^ W+W^ : iu+ d+) (^^1 + ^2 + ^3 l) 



J- 



01/ 2 \00 MOl 



(43) 



where (t± = {ai it i <T2)/\/2. Thus 



Atu. = + A3, AtcI. = ^Ai + A3. (44) 
We introduce now the isospin vector amplitude of Section [3] 

MA 

= M2 . (45) 




The Born amplitudes At1_ and a!j'^_ of Eq. ([36]) correspond to the vector 



^1 =2\l/u\. (46) 




For the amplitude A]^q_ of the left-handed quark-antiquark pair transition into Goldstone bosons, 
the isospin basis is: 

In the first term the sum over a goes from 1 to 3, but only (J3 contributes to the production of 
charged (j) pair so that 

Alu- = -jAa + At, A^d. = ^Aa + Ai. (48) 
As for the transverse case we introduce the isospin vector amplitude A-l of Section [3] 

A,. = (■^l) . (49) 
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The Born amplitudes Al^_ and A'^J_ of Eg. ([36]) correspond to the vector 



l(0) 



A! 



(0) 
■Lq- 



1 



^i^f 4.2 / .2 • 



2 / .,2\ ' 



12 



(50) 



where Yi (Yf) is the hypercharge of the initial (final) state. The isospin vector amplitudes can also 
be decomposed according to the gauge couplings 



A' 



(1) 

•P(jr_ 



A A! 



2 ^yl(l) 



■Pq-,SU(2) T V T V I "^-Pq-.QBB -^pg. 



(1) 



iP = T,L). (51) 



The expressions for the SU{2) and y components of the vectorial amplitude are then written in 
terms of the parameters of the evolution equation 



A' 



(1) 

Pq- ,i 



.(1) 



(i = 5C/(2),y). (52) 



The coefficients 7^^] , C^] and ^ are universal, the other are obtained by explicit calculation. 



/I) c^'^ 

lp,l ' Spjj — " ^p.i 

The expression for the amplitudes ^ and A^^^ ■ in terms of the anomalous dimensions can 
be then obtained by using Eqs. ipH liHj) and takes the form: 
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where we introduced the notations 
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The coefficient are given by: 
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Transverse W, left-handed quarks, Y component: 
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Longitudinal V^, left-handed quarks, SU(2) component: 
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Longitudinal W, left-handed quarks, Y component: 
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For the right-handed quarks the amphtudes are isospin singlet so that the SU (2) and Y contribu- 
tions factorize with respect to the Born ampUtude: 
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Longitudinal W^, right-handed quarks, SU(2) component: 
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Longitudinal W , right-handed quarks, Y component: 



1lu+,y 



A 



(1) 



A 



(1) 

Ld+,Y 



25 
18' 



1 4 



7 
3' 

209 25 

Ts" ~ 54' 



(1) 



X 



2 ^ 



4. 



.(1) . 



57 1 2Ar 4^2 



4 4 

13 
18' 



2 



c 



(1) 



Xz ■ 



_ 4 
94 13 



X 



(1) 



-L 



41 



— L_ LutL^ 1 vr ITT + Wy H — -L„ L 

6 3«t^ 9 54 6 2'' V6 4^ 4^/^ 9^ 



1 



1 



-T^ - - 

4 Xz 4 



(61) 



The NNLL two-loop amplitudes are obtained by using Eqs. (|9|44|48p : 

„(2) _ Irq- g(l) 2 
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where 

[PJ \pq^— g Jpq^,SU{2) g V-'fl?T>i'^ g ^y^fJq,QEDi J — 7)S,)A) [A^Ajpg-— g Anj., 51/(2)) [HX.\pq+—^^ 

= aSL,s.,), 411 = = 0, = 0; A^t = "i^* " 3^'' (63) 

with P = T,L and q = u,d. Note that the Yukawa contribution of Eq.psp is partiahy contained 
already in Cl ^su(2) remaining piece is given by ■ The square of the matrix x gives the 

fohowing contribution: 

r^2n(l) _ r 2l(l) , 2y^^^ ( Y^^^ - Y^^^ \ + ( Y^^^ - Y^^^ 

[A \Pq^ — [A \pq^L,SU(2) ^ ^Apq^,SU(2)\'^w M^QED I ^ 1 V^eT.^ ^,QED 

r^2n{l) _ r-2i(l) , ^Y^l) ( ,2 (l) _ „2 (1) A 

La \Pq^~ La \pq:p,SU(2) ' '^Apq^,SU(2)\''wAPqTf,Y '^y^r Xq,QED I 

(2) 

Finally, the pure two-loop quantities Tpq^ are given by: 
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(2) , 11 2 (2) _ 52 (2) , 9 2 (2) _ 520 
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lLq+,SU(2)— 3 "'"'^ ' 1lu+,Y — , ^Ld+,Y— ^-y ' 

Throughout this section the following notations have been used 
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Appendix B 

In this appendix we present the result for the one and two-loop corrections to the partonic cross 
section in NNLL approximation. The differential cross sections are obtained from the amplitudes 
given in the previous appendix through the relations: 
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The perturbative series Eq. (j24p for the cross section takes the form 
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We expand the corrections terms Jp*^^^ in powers of the large logarithm L = ln(s/M^) 
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Numerically for the one-loop coefficients we get 
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1.45- If - 0.45- li + 0.30 - 0.30 /„ + 1.01 ; 

u t 
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a 



(1) - _1 OK ^(1) - 



-1.25 , bl^^ = -0.43 It + 0.43 - 6.69 , 



c^l1+ = 0.20| ll - 0.20^ - 0.11 Z„ + 0.11 It + 46.23 ; 

a^^ = -1.20 , 6^+ = 0-22 k - 0.22 - 6.85 , 

It * t 



c['> = 0.10- 1'^ - 0.10- - 0.05 k + 0.05 Z„ + 46.31 . 



The two-loop coefficients read 

a^rL = 11-76 , b^rL = (^32.82 - 19.40^^ lu - 13.42 k - 17.34 , 



.(2) 



= 1^34.95 - 22.59^ + 16.04^ _ 10.72 Z| + 34.13 - 20.33^j/„Zt 



+ 1^- 44.44 + 37.21| J Z„ + 9.24 k - 28.46 + 39.48 

^{t — u)u 



19.56| + 12-04^ j ll + 14.55;^ I'f + ^- 22.06 + 24.08^^ /„ + 12.04 

aS5_ = 12.52, = (^36.94 - 20.01 16.93 /„- 20.89, 

c55_ = ^39.69 - 25.29^ + 16.42^^/^2 _ 9 28Z^ + ^- 41.00 - 18.48|^Z„Zt 



+ ( - 53.63 + 40. 13^ ) Zt + 15-58 Z„ - 23.56 + 39.48 * 

4 



20.18- + 12.42— Z; + 15.01^ It + - 22.77 + 24.85- k + 12.42 
u V? J ^ t \ u ' 



aiS_ = 3.12 , b^L = 12.42 - 2.42 k + 10.54 , 



^- 3.87^ + 12.34^ ll + (^1-37^ + 0.47^ Zf - 0.81 Z„ Z* + 23.93 Z„ - 6.34 Z* - 274.80 ; 



41 = 3.52 , 6f j_ = 13.76 Z* - 3. 14 + 9.44 , 

4d_ = (- 3-85^ + 13.44^/(2 + (^119^ + 0-70^/^ - 2.14 + 21.33 - 5.91 /„ - 101.71 ; 
ai'2+ = 0.78 , = 0.54 k - 0.54 /„ + 10.11 , 

cin+ = (oM + 0-25^ j ll + 0.16 - 0.25;^ j Zf - 0.19 lu k - 2.30 lu + 2.30 k - 88.00 ; 

«S+ = 0-72 , = -0.26 k + 0.26 + 10.03 , 

= (0.14 + 0.12|^Zf + ^- 0.10 - 0.12^^/2 _ 0.05 - 1.19 Zt + 1.19 Z„ - 85.17. 

Here lu = ln{—u/s) and It = ln(— t/s). 
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Appendix C 

In Ref. ^21j the contribution of the imaginary part of the anomalous dimension matrix x^^^ 
(given in Eq. ()55p above) has been missed in the numerical estimates. This contribution changes 
the NNLL two-loop correction in the transverse boson production cross section. It results in an 
additional term 

in the coefficient of the quadratic logarithm in Eqs. (14, 33) of Ref. |21] . 
References 

[1] V.S. Fadin, L.N. Lipatov, A.D. Martin, and M. Melles, Phys. Rev. D 61 (2000) 094002. 

[2] J.H. Kiihn, A.A. Renin, and V.A. Smirnov, Eur. Phys. J. C 17 (2000) 97; Nucl. Phys. B (Proc. 
Suppl.) 89 (2000) 94. 

[3] M. Melles, Phys. Rev. D 63 (2001) 034003; D 64 (2001) 014011. 

[4] A. Denner, M. Mehes, and S. Pozzorini, Nucl. Phys. B 662 (2003) 299. 

[5] J.H. Kiihn, S. Moch, A.A. Renin, and V.A. Smirnov, Nucl.Phys. B 616 (2001) 286, Erratum 
ibid. B 648 (2003) 455. 

[6] B. Feucht, J.H. Kiihn, A.A. Renin, and V.A. Smirnov, Phys. Rev. Lett. 93 (2004) 101802. 

[7] B. Jantzen, J.H. Kiihn, A.A. Renin, and V.A. Smirnov, Phys. Rev. D 72 (2005) 051301 (R); 
Nucl. Phys. B 731 (2005) 188, [Erratum ibid. B 752 (2006) 32. 

[8] J. H. Kiihn and A. A. Renin, |arXiv:hep-ph/9906545[ 

[9] M. Ciafaloni, P. Ciafaloni and D. Comehi, Phys. Rev. Lett. 84 (2000) 4810. 

[10] A. Denner and S. Pozzorini, Eur. Phys. J. C 18 (2001) 461; A. Denner and S. Pozzorini, Eur. 
Phys. J. C 21 (2001) 63; S. Pozzorini, Nucl. Phys. B 692 (2004) 135. 

[11] J. Chiu, A. Fuhrer, R. Kelley, A.V. Manohar, Phys. Rev. D 77 (2008) 053004. 

[12] M. Lemoine and M.J.G. Veltman, Nucl. Phys. B 164 (1980) 445. 

[13] M.Bohm, A.Denner, T.Sack, W.Beenakker, F.A. Berends, and H. Kuijf, Nucl. Phys. B 304 
(1988) 463. 

[14] J. Fleischer, F. Jegerlehner, and M. Zralek, Z. Phys. C 42 (1989) 409. 

[15] W. Beenakker, A. Denner, S. Dittmaier, R. Mertig and T. Sack, Nucl. Phys. B 410 (1993) 
245. 

[16] W. Beenakker, F.A. Berends, A.R Chapovsky, Nucl. Phys. B 548 (1999) 3. 

[17] S. Jadach et al., Comput. Phys. Commun. 140 (2001) 432; Phys. Rev. D 65 (2002) 093010; 



23 



[18] A. Denner, S. Dittmaier, M. Roth and D. Wackeroth, Nucl. Phys. B 587 (2000) 67; Comput. 
Phys. Commun. 153 (2003) 462. 

[19] A. Denner, S. Dittmaier, M. Roth and L.H. Wieders, Phys. Lett. B 612 (2005) 223; Nucl. 
Phys. B 724 (2005) 247. 

[20] M. Beccaria, F.M. Renard, and C. Verzegnassi, Nucl. Phys. B 663 (2003) 394. 

[21] J. H. Kiihn, F. Metzler and A. A. Penin, Nucl. Phys. B 795 (2008) 277 [Erratum-ibid. 818 
(2009) 135]. 

[22] E. Accomando, A. Denner, and S. Pozzorini, Phys. Rev. D 65 (2002) 073003. 

[23] E. Accomando, A. Denner, and A. Kaiser, Nucl. Phys. B 706 (2005) 325. 

[24] J. Chiu, A. Fuhrer, R. Kelley, A.V. Manohar, Phys. Rev. D 80 (2009) 094013. 

[25] U. Baur, T. Han, and J. Ohnemus, Phys. Rev. D53 (1996) 1098. 

[26] L.J. Dixon, Z. Kunszt, and A. Signer, Phys. Rev. D 60 (1999) 114037. 

[27] M. Grazzini, JHEP 0601 (2006) 095. 

[28] T. Binoth, M. Ciccolini, N. Kauer et al, JHEP 0612 (2006) 046. |hep-ph/0611170] . 

[29] V.V. Sudakov, Zh. Eksp. Teor. Fiz. 30 (1956) 87. 

[30] R. Jackiw, Ann. Phys. 48 (1968) 292; 51 (1969) 575. 

[31] A.H. Mueller Phys. Rev. D 20 (19) 2037 

[32] J.C. Collins, Phys. Rev. D 22 (1980) 1478; Adv. Ser. Direct. High Energy Phys. 5 (1989) 573. 

[33] A. Sen, Phys. Rev. D 24 (1981) 3281. 

[34] A. Denner, B. Jantzen and S. Pozzorini, Nucl. Phys. B 761 (2007) 1. 

[35] J. Frenkel and J.C. Taylor, Nucl. Phys. B 116 (1976) 185. 

[36] A. Sen, Phys. Rev. D 28 (1983) 860. 

[37] G. Sterman, Nucl. Phys. B 281 (1987) 310. 

[38] J. Botts and G. Sterman, Nucl. Phys. B 325 (1989) 62. 

[39] M. Melles, Phys. Rev. D 64 (2001) 014011. 

[40] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Phys. Lett. B 531 (2002) 216. 

[41] T. Hahn, Comput. Phys. Commun. 168 (2005) 78. 

[42] M.Bohm, A. Denner and H. Joos, "Gauge theories of the strong and electroweak interaction", 
Teubner, Stuttgart (2001). 



24 



